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A shortcut to adiabaticity is a finite-time process that produces the same final state as would result 
from infinitely slow driving. We show that such shortcuts can be found for weak perturbations from 
linear response theory. With the help of phenomenological response functions, a simple expression 
for the excess work is found-quantifying the nonequilibrium excitations. For two specific examples, 
i.e., the quantum parametric oscillator and the spin 1/2 in a time-dependent magnetic field, we show 
that finite-time zeros of the excess work indicate the existence of shortcuts. Finally, we propose a 
degenerate family of protocols, which facilitates shortcuts to adiabaticity for specific and very short 
driving times. 
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I. INTRODUCTION 

Thermodynamics is a phenomenological theory to de¬ 
scribe the transformation of heat into work. How¬ 
ever, only quasistatic, i.e., infinitely slow processes are 
fully describable by means of conventional thermody¬ 
namics [1]. For all realistic, fmite-tim.e-nonequilibrium- 
processes, the second law of thermodynamics constitutes 
merely an inequality, expressing that some portion of 
the energy or entropy is irreversibly lost into nonequi¬ 
librium excitations. For isothermal processes, this “loss” 
is quantified by the excess work ( W ex ), which is the differ¬ 
ence between the total nonequilibrium work ( W) and the 
work performed during a quasistatic-egm/z&rmm-process 
(Wq S ), (Wex) = (W) — (Wq S ). For macroscopic, open sys¬ 
tems, (Wq S ) is simply given by the free energy difference 
A F. However, the identification of the equilibrium work, 
(Wq S ), with the free energy difference, A F, is only true 
for open systems. For isolated systems, the minimal work 
is not given by the free energy difference and (W qs ) has to 
be analyzed carefully [2]. In addition, for quantum sys¬ 
tems, the situation is particularly involved as quantum 
work is not an observable in the usual sense, as there is 
no Hermitian operator, whose eigenvalues are given by 
the classical work values [3-7]. 

Nevertheless, finding “optimal” quantum processes, 
for which only the minimal amount of (W e x) is lost 
into nonequilibrium excitations, is of fundamental im¬ 
portance. Consequently, a lot of theoretical and experi¬ 
mental research has been dedicated to the design of so- 
called shortcuts to adiabaticity, i.e., finite-time processes 
with suppressed nonequilibrium excitations [8]. To this 
end, a variety of techniques has been proposed: the use 
of dynamical invariants [9], the inversion of scaling laws 
[10], the fast-forward technique [11, 12], and transition¬ 
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FIG. 1. (Color online) Sketch of the thermodynamic processes 
under study. At t = to, the system is prepared in equilibrium 
with inverse temperature /3 before the system is decoupled 
from the environment and controlled externally from t = to + 
0+ until a final time tf. 


less quantum driving [13-16]. All methods have in com¬ 
mon that practical implementations are rather involved 
as the full dynamics has to be solved to determine the 
shortcut. Therefore, more recent research efforts have 
been focusing on identifying optimal protocols from op¬ 
timal control theory [17, 18], properties of the quantum 
work statistics [19], or “environment” assisted methods 
[ 20 ]. 

The present analysis is dedicated to finding shortcuts 
to adiabaticity from a phenomenological approach-linear 
response theory. For classical systems, it has been re¬ 
cently shown that there exist finite-time processes with 
zero excess work [21], In this paradigm, (W ex ) is fully 
determined by the phenomenological response of the sys¬ 
tem to an external perturbation [22, 23]. Thus, we nei¬ 
ther have to solve the dynamics [13-16] nor do we have 
to determine the quantum work statistics [19] to mini¬ 
mize (W ex ) • In the following, we will extend our previous 
findings [21, 23] to the quantum domain. To this end, we 
will consider a thermally isolated quantum system under 
weak perturbation and derive a linear response expres- 









2 


sion for (W e x)- After establishing the general theory, we 
will turn to analytically solvable and pedagogically elu¬ 
cidating examples, namely the parametric harmonic os¬ 
cillator and the spin 1/2 in a time-dependent magnetic 
field. This will allow us to study the range of validity 
of the linear response approach by comparing our find¬ 
ings with the exact results from the full quantum work 
statistics [24, 25]. We will show that the protocols with 
zero excess work from linear response theory, indeed, they 
facilitate transitionless quantum driving for weak pertur¬ 
bations. Finally, we will propose a family of degenerate 
protocols, which facilitates shortcuts to adiabaticity for 
arbitrarily fast driving. 

II. QUANTUM WORK FROM LINEAR 

RESPONSE THEORY 

We begin by generalizing the previous classical treat¬ 
ment of the excess work (W ex ) [21] to the quantum do¬ 
main. Imagine a quantum system with time-dependent 
Hamiltonian H t , which is prepared initially in a thermal 
equilibrium state, po = exp(—/3Ho)/Zo, where Zq is the 
partition function, Z 0 = tr{exp (— (3H a )}. At t = fo+0 + , 
the system is decoupled from the environment and the 
Hamiltonian is varied according to some protocol A t with 
H t = H( A t ). Such a processes is sketched in Fig. 1. 

The external control parameter At is written as 

At = Ao + SX g(t), (1) 

where At starts in an initial value Ao, <5 A is the amplitude, 
and g(t) obeys g(t 0 ) = 0 and g(tf) = 1. Thus, A t varies 
from Ao to A/ = Ao + <5A. 

For small systems, work is a fluctuating quantity [26] 
and for a specific protocol g(t) the average work reads 

(W) = f ' diAt (d x H) , (2) 

Jto 

where the angular brackets denote an average over many 
realizations of the same process and the dot denotes a 
derivative with respect to time. 

We will now evaluate the general expression for the 
average work (2) by means of linear response theory. To 
this end, we expand the Hamiltonian up to linear order 
in the amplitude <5A, 

H{ A t ) = H( Ao) + <5A g(t) d x H + 0(6 A 2 ). (3) 

By substituting Eq. (3) into Eq. (2) and identifying d\H 
as the generalized force [21, 23, 27, 28], it can be shown 
[21] that the average work (2) becomes 

(W) = s\ (dxH) po ~ { -^*( 0 ) 

rtf pt—to ( 4 ) 

- {SX) 2 dtd t g d s^f(s)d a g(t - s), 

J to J 0 

where 'F(t) is the relaxation function [27, 28]. 


Until Eq. (4), the present treatment is identical to the 
classical case [21]. However, in the quantum case, the 
relaxation function ’F(t) is determined by the quantum 
response function <j>(t), 4>(t) = —’F(t), with [27, 28] 

<t>(t) = tr {p 0 [A 0 , A t ]} , (5) 

where A = d\H is the generalized force. To avoid clutter 
in the formulas, we introduce in Eq. (5) the notation 
A(t) = At- 

In complete analogy to the classical case [21], the first 
two terms of Eq. (4) are independent of the specific proto¬ 
col g(t ) and we identify the quasistatic, equilibrium work 
as 

<Wqs> = 5X ( d x H) po - M!*(0). (6) 

In the remainder of this analysis, we will analyze the 
excess work, 

(We X ) =-(<5A) 2 / dtd t g [ ds^(s)d s g(t-s) (7) 
J to J 0 

for two analytically solvable examples. We will show that 
whenever this thermodynamic quantity vanishes in finite 
time, the quantum adiabatic invariant is conserved and 
therefore the system can be driven through a shortcut to 
adiabaticity. 

Generally, it is easy to see that if the adiabatic theorem 
is fulfilled, no transitions between eigenstates occur, and 
therefore the excess work (W e x) has to vanish. However, 
the reverse is not necessarily true. Even if the excess work 
vanishes, one could imagine a process during which some 
transitions between eigenstates do occur, however in such 
a way that their energetic contribution “cancels out.” In 
the following, we will analyze this issue with the help of 
two fully analytically solvable examples-the parametric 
harmonic oscillator and a spin-1/2 particle in a magnetic 
field. We will find that at least within the range of valid¬ 
ity of linear response theory, such “canceling” transitions 
do not occur since for a “shortcut” not only the excess 
work vanishes, but also the adiabatic invariant is (ap¬ 
proximately) conserved. For classical systems, a similar 
analysis was developed in Ref. [21]. 

III. PARAMETRIC HARMONIC OSCILLATOR 

We consider the time-dependent Hamiltonian 

H (X t ) = \ + 2 Xtx2 

where x and p are the coordinate and momentum opera¬ 
tors, respectively. This system can be solved analytically 
[24, 25] for specific protocols A t that drive the system 
from an initial to final value of A, as illustrated in Fig. 2. 
To simplify notation, we further set to = 0 and tf = r. 




3 



x 


B. Exact solution 

The parametric harmonic oscillator (8) has been ex¬ 
tensively studied, since it can be solved analytically 
[24, 25, 29, 30] for specific driving protocols and it de¬ 
scribes quantum thermodynamic experiments in cold ion 
traps [31-33]. The time-dependent mean energy can be 
written as [25, 34] 


(Hr) 


2 


Q* coth 




(13) 


where <5* is a measure of adiabaticity [24, 25, 29]. This 
measure is fully determined by two special solutions, X t 
and Y t , of the force-free equation of motion [29], 


FIG. 2. (Color online) Parametric harmonic oscillator (8) 
with Ao (dashed line) at time t = to and A/ (solid line) at 
t — T. 


A. Linear response approach 

The response function (5) is obtained by solving 
Heisenberg’s equations of motion for a fixed, initial value 
of A. Hence, we obtain, after a few simple lines, 

(f){t) = coth sin(2-\/Ao t). (9) 

It is interesting to note that the system’s response is oscil¬ 
latory. Consequently, we have the “relaxation” function 


x t + X t x t = 0. 


(14) 


We have 


Q* 


1 

2a/AoA / 


Ao (A f X 2 + X 2 ) + (Ay Y 2 + Y 2 )] , 

(15) 


with X 0 = 0, X 0 = 1 and Y 0 = 1, Y 0 = 0. [29]. Note that 
these initial conditions for X t and Y t are chosen for the 
sole sake of simplifying the mathematical treatment [29] . 
For the quantum harmonic oscillator, the time-dependent 
action S = E(t)/u(t) is conserved if [25] 


X 2 + X t X? 
\/A t 



Y t 2 + ^tY t 2 
\/A t 


\/Ao ■ (16) 


= 2\^/x q COth ( /3ri ^ A ° ) cos (2a/Ao i) • (10) 

Generally, relaxation functions describe how a system re¬ 
laxes towards an equilibrium state. However, since the 
present system has only a single degree of freedom and 
is thermally isolated, the “relaxation” function exhibits 
nondecreasing oscillations. 

For the sake of simplicity, we further assume that the 
stiffness varies linearly with time, 

A* = Aq + SX t/r , (11) 


for which we obtain 



( 12 ) 

Equation (12) constitutes our first main result. In com¬ 
plete analogy to the classical case [21], the excess work 
vanishes for all zeros of the sine function, i.e., for all 
r = nn/y/Xo with n being an integer. In the classical 
case, these “special” driving times have been attributed 
to a conservation of the adiabatic invariant during the 
finite-time process [21]. 

In the next section, we will further analyze this obser¬ 
vation and show that the minima of (W ex ) (12), indeed, 
identify shortcuts to adiabaticity. 


Thus, it is easy to see that Q* > 1, where the equality 
holds for quasistatic processes. Accordingly, the exact 
expression for the excess work reads 

(w e r ct ) = RR±R co th ( q* -1). 

(17) 

Note that Q* depends only implicitly on the protocol X t 
through the solutions of Eq. (14). Therefore, it is ad hoc 
not clear whether the exact excess work (17) exhibits the 
same zeros as the expression from linear response theory 
(12) for the linear protocol (11). 

To gain insight and to build intuition, we plot the mea¬ 
sure of adiabaticity Q* (15) for the linear protocol in 
Fig. 3 for various strengths of the perturbation 5X. We 
observe that generally Q* — 1 exhibits oscillations, but 
no zeros as a function of r. 

For weak driving, however, i.e., 5X/Xq <C 1, where we 
expect linear response theory to hold, the minima of Q* — 
1 get infinitely close to zero. In Fig. 4, we compare the 
excess work from linear response theory (12) with the 
behavior of Q* — 1 for weak driving. We observe very 
good agreement between the result from linear response 
theory (12) and Q* — 1. 

It has also been shown that for Q* = 1, the quan¬ 
tum adiabatic theorem is fulfilled, i.e., for such processes 
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FIG. 3. (Color online) Measure of adiabaticity Q* (17) as a 
function of the switching time for the linear protocol, g(t) = 
t/r, and A/ = 2.0 (blue solid line), A/ = 1.7 (green dashed 
line), Xf = 1.5 (yellow dot-dashed line), A/ = 1.3 (orange 
dash-dotted line), and A/ = 1.1 (red dotted line), and Ao = 
1 . 0 . 



FIG. 4. (Color online) Excess work from linear response the¬ 
ory (12) (blue solid line) together with Q*- 1 (black dashed 
line) as a function of r for the linear protocol (11) and 
<5 A = 0.1. The symbol W EX denotes {W EX ) measured in units 
of (<5A/^A4) 2 {hV AO) coth (/3hVT q/2) /4. 


there are no transitions between different energy eigen¬ 
states [29]. Thus, we conclude that the zeros of the ex¬ 
cess work, indeed, identify finite driving times for which 
transitionless quantum driving is facilitated-shortcuts to 
adiabaticity from linear response theory. 


C. Range of validity of linear response theory 

Linear response theory can be understood as a phe¬ 
nomenological theory of weak perturbations [27]. Thus, 


the numerical and qualitative agreement between exact 
(17) and approximate (12) results cannot be considered 
satisfactory. To deepen the insight into the approxima¬ 
tions, we will now derive Eq. (12) from the exact expres¬ 
sion (17) without having to rely on phenomenology. 

To this end, we expand the exact expression (17) in 
powers of <5A up to second order. Note that Q* depends 
implicitly on the protocol X t and we write Q*((5A). We 
have 

(^exr) * \ COth jjA yfcd X Q*(0) 

+ [c> a Q*(0) + Ao^Q*( 0)] +o(5A 3 )j, 

(18) 

where we used Q*(0) = 1. We now have to show that 
there exist approximate solutions X t and 34 of the equa¬ 
tion of motion (14) such that Eq. (18) reduces to the 
linear response expression (12) with X t and 34 replacing 
X t and Y t in Eq. (15). 

Comparing Eqs. (12) and (18), we conclude that X t 
and 34 have to fulfill 

d\Q*(0) = 0 and d 2 Q*(0) = Sm ^ 2 ° T) . (19) 

Additionally, we know that X t and 34 have to obey 
X t y t — X t y t = 1 [29]. The latter condition is just an 
expression of the commutation relation between position 
and momentum [29]. For S A = 0, the solution of Eq. (14) 
is given by the sine and cosine function [29]. Hence, we 
make the ansatz 

X t = — L= sin + SX F t + 0(SX 2 ) 

VA 0 v > (20) 

34 = cos (^\Z~Xf)t^ + SXQt + C3(4A 2 ), 


where T t and Qt are two time-dependent functions deter¬ 
mined by the conditions (19). 

It is then a tedious but straightforward exercise to show 
that 


Tt = 


t 2 + 4aAoT 
4A 0 r 


cos 



and 


t — 4bX 0 r 

-sin 

4Aq vA qt 



Qt 


t 2 + 4aAoT 
4A 0 r 


cos 



+ 


t 2 Ao — 4 cAot — 1 
4Xo^/Xqt 


sin 



( 22 ) 


The three constants a, b , and c are determined by the 
boundary conditions JF 0 = a, Fq = b and Qq = — 6 , Qo = c 
[35]. The expressions of T and Q are rather lengthy and 
can be found in Appendix A. 

The solutions (20) together with Eqs. (21) and (22) 
are the approximate solutions of Eq. (14), for which the 
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FIG. 5. (Color online) (a) Exact solution Xt (solid line) and 
approximate solution Xt (20) for SX = 0.1 (dashed line), (b) 
Exact solution Y t (solid line) and approximate solution y t 
(20) for S X = 0.1 (dashed line). Shaded area signifies a SX 
environment around the exact results. 


exact expression for the excess work (17) reduces to the 
result from linear response theory (12). In Fig. 5, we plot 
the approximate solutions (20) together with the exact 
solutions of (14). We observe that X t and W are within 
a 5 A environment around the exact results, as one would 
intuitively expect by construction. 

In conclusion, we have shown that results from lin¬ 
ear response theory can also be obtained from expanding 
the exact solutions for weak driving. Thus, the linear 
response expressions are not only considered to be quali¬ 
tatively and phenomenologically true, but also quantita¬ 
tively exact. 


D. Optimal protocols — shortcuts to adiabaticity 


In an analogous classical treatment, it has been shown 
that the linear parametrization (11) is not the only proto¬ 
col with zero excess work. Rather, there is a degenerate 
family of optimal protocols [21, 36] for which nonequi- 
librium excitations are suppressed. This family is given 



\Ao r 


FIG. 6. (Color online) Excess work (12) (black dashed 
line) and normalized adiabatic parameter Q* — 1 (red solid 
line) as a function of the switching time for the optimal 
protocols (23) with (a) a = 1, k = 2 and (b) a = 1, 
k = 4. The symbol W ex denotes (W ex ) measured in units 
of ((5A/x/A[[) 2 (hVX,) coth {phV\o/2) /4. 


by 


g{t) = t/r + a sin(K7rf/r) , (23) 


where k is an integer and a is any arbitrary real number. 

The quantum excess work (12) merely differs in the 
prefactor from the classical expression 


<Wex) 


npV a 7 «i 


SX \ 2 1 sin 2 (v / Ao r) 
2/3 Aq t 2 


(24) 


which is obtained in the limit H/3\/Xo <C 1. Thus, the 
degenerate class (23) constitutes a family of shortcuts to 
adiabaticity for the quantum harmonic oscillator under 
weak driving. Figure 6 illustrates (W ex ) (12) together 
with Q* — 1 for two members of the family (23). It has 
been shown [21] that the shortcut to adiabaticity is ob¬ 
tained for v / Ao r = wr, with n integer, and 


V^o T 


(ktt/2) 

(1 + Kna) 1 / 2 


(25) 
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Finally, it is worth emphasizing that such shortcuts to 
adiabaticity can be obtained for arbitrarily short switch¬ 
ing times by choosing a appropriately [21]. 

IV. SPIN 1/2 IN A TIME-DEPENDENT 
MAGNETIC FIELD 

Our second example is a spin 1 /2 in a time-dependent 
magnetic field subjected to the constraint |B(t) | = B 0 = 
constant. Its Hamiltonian reads 

= fc<r- B( i), (26) 

where er denotes the Pauli matrices. Due to the above- 
mentioned constraint on B(i), it is more convenient to 
choose the following parametrization 

/sin [</>(£)] cos [d(f)]\ 

B(t) = B 0 sin [^(t)] sin [9(t.)} . (27) 

V cos [<p{t)} ) 

Hence, the time dependence of the set of allowed pro¬ 
cesses parameterized by the angles ip(t) and 0{t) is, in 
analogy to Eq. (1), expressed as 

ip(t) = ip 0 + 5ipg v (t), (28a) 

0(t) = 9 0 + 59 ge(t) (28b) 

where the boundary conditions g v ,e(0) = 0 and g Vl e(r) = 
1 must hold. 

Linear response theory provides a good description of 
( W ex ) as long as 6ip and 69 are sufficiently small. In 
this regime, one can easily show that the angle 9(t) plays 
no role and the thermodynamic work (7) depends on the 
nonequilibrium of d v H only. Thus, the response function 
is given by Eq. (5) with A t = d^H(t) and it is straight¬ 
forward to obtain 

\ (7-®o) 2 tanh sin (jB 0 t) , (29) 

from which, using again <p(t) = — 'L(f), we have the re¬ 
laxation function 

^ tanh ^ —cos (u Q t) , (30) 




CUqT 


FIG. 7. Time evolution of the real and imaginary parts of the 
coefficients c+(t) and c_(t) given by Eq. (32) for the initial 
condition c+(0) = 1 and c_(t) for cos (5ipt/2r) ~ 1. 


A. Quantum adiabatic invariant 

Analogously to Sec. Ill, we will now verify that the 
quantum adiabatic invariant is conserved for spin-1/2 
particles driven by Eq. (23). To this end, we analyze 
the time evolution of the coefficients c_|_ (t) and c_ (t) ap¬ 
pearing in the expansion 


where we defined ujq = 7-Bo- 

The time dependence of the relaxation functions (10) 
and (30) has the same functional form. Therefore, the ex¬ 
cess work performed by an external agent while driving 
the spin 1/2 will behave exactly the same as in the para¬ 
metric harmonic oscillator. For instance, the protocols 
given by (23) also constitute a family of optimal proto¬ 
cols for the present system. Nevertheless, the values of 
ojqt for which the excess work vanishes are a bit different 
from those in Fig. 4 due to the absence of the factor 2 in 
cos (w 0 t) of Eq. (30). The linear protocol generates zeros 
for uiqt = n 2n. 


^ c„(f)exp 

n=+,- 



d t'E n (t'] 


\n;t ), 


(31) 

of an arbitrary state \ip(t)). We denote by E n (t ) and | n; t) 
the instantaneous eigenvalues and eigenstates of (26). 
The quantum adiabatic invariant is then conserved in fi¬ 
nite time if, after starting with c+(0) = 1 and c_(0) = 0 
at the beginning of a certain protocol g v (t), we obtain 
c+(t) = c+(0) and c_(r) = c_(0). 

The equations of motion for c+,_(f) are easily de¬ 
rived following standard procedures [37-39]. For the 
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parametrization (27) of B(f), we obtain 

dC ^ = ~^ T cos exp ^ iuj ot) c-(t), (32a) 

dC ~ d ^ = cos exp MO ’ ( 32b ) 

considering $o = 0, ge{t) = 0, and g v {t) = t/r. 

Figure 7 shows the real and imaginary parts of the 
solutions of (32) as functions of ojqt for cos (6ip t/2r) ~ 
1, since we are in the regime Sip <Cl (see Appendix B 
for their analytical form). Since the initial conditions 
are c + (0) = 1 and c_(0) = 0, we should have a finite¬ 
time conservation of the adiabatic invariant every time 
we get a recurrence of this values. In Fig. 7, we see 
that this holds true for uj 0 t = n 2tt, although due to 
our approximations the imaginary part of c_ (t) does not 
vanish at these values of r. 


V. COMPLEX SYSTEMS 

The two case studies in Secs. Ill and IV are analytically 
solvable and pedagogically elucidating. In particular, we 
obtained exact expressions for the response functions (9) 
and (29). However, this is not feasible for general and 
more realistic systems with more degrees of freedom. It 
has been shown [22, 23, 27, 28] that linear response theory 
performs well when only phenomenological information is 
known about the system of interest. In other words, even 
when the response function is not exact, the predictions 
of linear response theory provide good approximations. 
Finding shortcuts from linear response theory and opti¬ 
mizing {W ex ) circumvents the difficult problem of having 
to solve for the full quantum dynamics. 

It has been shown that ( W ex ) will have finite-time min¬ 
ima, or nonmonotonic behavior as a function of r, if the 
relaxation function is sufficiently oscillatory. This can be 
illustrated, for instance, using the following phenomeno¬ 
logical ansatz [23, 28]: 

'F(f) = 'F(O) exp(—a|t|) ^cos(wt) + — sin (cut) ^ , (33) 

for the relaxation function. Plugging the expression 
above in Eq. (7), we obtain the results shown in Fig. 8 for 
different values of a/uj. As this ratio decreases, the excess 
work starts to show minima whose value approaches zero. 
There are several systems for which Eq. (33) describes the 
relaxation dynamics very well. Among them, we mention 
a system composed of weakly interacting magnetic mo¬ 
ments in the regime where Bloch equations are valid [40] . 


VI. CONCLUDING REMARKS 


Identifying optimal quantum processes with sup¬ 
pressed or even vanishing nonequilibrium excitations is 



0 2 4 6 8 10 12 


CUT 

FIG. 8. (Color online) Excess work Wf/ in units of 
(5ip) 2 'S>(0)/2 for the linear protocol g v (t) = t/r and the re¬ 
laxation function (33). 


an important topic, which has recently been attracting 
intense research efforts. However, all methods currently 
available necessitate the solution of the full quantum dy¬ 
namics. In the present work, we have proposed a phe¬ 
nomenological alternative. By generalizing our previous 
result for the excess work from linear response theory 
to a quantum system, we have shown that shortcuts 
to adiabaticity can be identified from a mathematically 
simple theory. This observation has been proven for 
two paradigmatic examples of quantum thermodynam¬ 
ics, namely the parametric harmonic oscillator and the 
spin 1 /2 in a time-dependent magnetic field. 
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Appendix A: Approximate solution within linear response 


The full expressions for the approximate solutions X t and 34 in Eqs. (21) and (22) are given in terms of the three 
constants a, b , and c. These are determined by solving the force-free equation of motion (14) with the boundary 
conditions Tq = a, Tq = b and Q o = —Go = c. We have 


— 1 — 2A 0 t 1 2 3 4 5 6 7 8 + cos(2v4\ot) — 2-/Ao'rsin(2v / AoT) _ ^ —1 + 2A 0 t 2 + cos(2v4\ot) — 2^/X^t sin(2\/Aor) 


8A§r 


8A 0 r 


and 


b = 


8Ao-\/Ao t 


2 + \J 4 + 2A 0 r 2 - 4 cos(2a/A^t) + 2a/Aqt cos(2a/Aqt) + sin(2i/AoT) 


(Al) 

(A2) 


Appendix B: Derivation of time-dependent coefficients c+ and c_ 

According to Ref. [39], the coefficients c + and c_ satisfy the differential equations for Eq. (27), 

dC+ ^ = _ 47 cos (jfr*) exp (-iw 0 i)c_(i), 

dc-(t) Sp f Sp \ , . 

^ c °s ( ) exp (*w 0 f)c + (f). 


(Bla) 

(Bib) 


In the regime <C 1, we make the approximation cos[J(^/2rt] ~ 1. Next, we solve exactly the equations with the 
initial conditions c+(0) = 1 and c_(0) = 0. After making f = r, we obtain the following equations: 


c + (w 0 t) = 


e » 


Sp 2 + 4 (uj 0 t) 2 


(Sp 2 + 4(w 0 r) 2 ) cosh ( - \J~(8p 2 +4(w 0 r) 2 ) 


- 2iu>oT\J —(Sp 2 + 4(w 0 t) 2 ) sinh ( - sj-(Sp 2 + 4(w 0 t) 2 ) 


(B2) 


c_(w 0 r)= 


e 2 




V + iM* “ nh ( 3 V-(V + 4 


where ojq = 'yB. 


(B3) 


[1] H. B. Callen, Thermodynamics and an Introduction to 
Thermostatistics (Wiley, New York, 1985). 

[2] A. E. Allahverdyan and T. M. Nieuwenhuizen, Phys. Rev. 
E 71, 046107 (2005). 

[3] H. Tasaki, arXiv:cond-mat/0009244. 

[4] J. Kurchan, arXiv:cond-mat/0007360. 

[5] P. Talkner, E. Lutz, and P. Hanggi, Phys. Rev. E 75, 
050102 (2007). 

[6] M. Campisi, P. Hanggi, and P. Talkner, Rev. Mod. Phys. 
83, 771 (2011). 

[7] P. Hanggi and P. Talkner, Nat. Phys. 11, 108 (2015). 

[8] E. Torrontegui, S. S. Ibanez, S. Martinez-Garaot, 
M. Modugno, A. del Campo, D. Guery-Odelin, 
A. Ruschhaupt, X. Chen, and J. G. Muga, Adv. At. 
Mol. Opt. Phys. 62, 117 (2013). 


[9] X. Chen, A. Ruschhaupt, S. Schmidt, A. del Campo, 
D. Guery-Odelin, and J. G. Muga, Phys. Rev. Lett. 104, 
063002 (2010). 

[10] A. del Campo and M. G. Boshier, Sci. Rep. 2, 648 (2012). 

[11] S. Masuda and K. Nakamura, Phys. Rev. A 84, 043434 

( 2011 ). 

[12] E. Torrontegui, S. Martmez-Garaot, A. Ruschhaupt, and 
J. G. Muga, Phys. Rev. A 86, 013601 (2012). 

[13] M. Demirplak and S. A. Rice, J. Chem. Phys. A 107, 
9937 (2003). 

[14] M. Demirplak and S. A. Rice, J. Chem. Phys. B 109, 
6838 (2005). 

[15] M. Berry, J. Phys. A: Math. Theor. 42, 365303 (2009). 
1161 S. Deffner, C. Jarzynski, and A. del Campo, Phys. Rev. 

X 4, 021013 (2014). 

[17] D. Stefanatos, Automatica 49, 3079 (2013). 
























9 


[18] S. Campbell, G. De Chiara, M. Paternostro, G. M. 
Palma, and R. Fazio, Phys. Rev. Lett. 114, 177206 
(2015). 

[19] G. Xiao and J. Gong, Phys. Rev. E 90, 052132 (2014). 

[20] S. Masuda and S. A. Rice, arXiv:1410.1883. 

[21] T. V. Acconcia and M. V. S. Bonanga, Phys. Rev. E 91, 
042141 (2015). 

[22] D. A. Sivak and G. E. Crooks, Phys. Rev. Lett. 108, 
190602 (2012). 

[23] M. V. S. Bonanga and S. Deffner, J. Chem. Phys. 140, 
244119 (2014). 

[24] S. Deffner and E. Lutz, Phys. Rev. E 77, 021128 (2008). 

[25] S. Deffner, O. Abah, and E. Lutz, Chem. Phys. 375, 200 

( 2010 ). 

[26] C. Jarzynski, Nat. Phys. 11, 105 (2015). 

[27] R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

[28] R. Kubo, M. Toda, and N. Hashitsume, Statisti¬ 
cal Physics II - Nonequilibrium Statistical Mechanics 
(Springer-Verlag, Berlin, 1985). 

[29] K. Husimi, Prog. Theor. Phys. 9, 381 (1953). 


[30] Z. Gong, S. Deffner, and H. T. Quan, Phys. Rev. E 90, 
062121 (2014). 

[31] G. Huber, F. Schmidt-Kaler, S. Deffner, and E. Lutz, 
Phys. Rev. Lett. 101, 070403 (2008). 

[32] O. Abah, J. Rofinagel, G. Jacob, S. Deffner, F. Schmidt- 
Kaler, K. Singer, and E. Lutz, Phys. Rev. Lett. 109, 
203006 (2012). 

[33] J. Rofinagel, O. Abah, F. Schmidt-Kaler, K. Singer, and 
E. Lutz, Phys. Rev. Lett. 112, 030602 (2014). 

[34] S. Deffner and E. Lutz, Phys. Rev. E 87, 022143 (2013). 

[35] Note that the initial conditions for F and Q are an ansatz 
chosen to fulfill Eq. (19). 

[36] T. Schmiedl, E. Dieterich, P.-S. Dieterich, and U. Seifert, 
J. Stat. Mech.: Theor. Expt. 2009, P07013 (2009). 

[37] M. V. Berry, Proc. R. Soc. London A 392, 45 (1984). 

[38] A. Bohm, Quantum Mechanics: Foundations and Appli¬ 
cations (Springer-Verlag, New York, 1994). 

[39] A. C. Aguiar Pinto, M. C. Nemes, J. G. Peixoto de Faria, 
and M. T. Thomaz, Am. J. Phys. 68, 955 (2000). 

[40] R. M. White, Quantum Theory of Magnetism (Springer, 
Berlin, 2006). 


